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We investigate the behavior of large amplitude oscillations 
of a trapped Bose-condensed gas of alkali atoms at zero tem- 
perature, by solving the equations of hydrodynamics for col- 
lective modes. Due to the atom-atom interaction, the equa- 
tions of motion are nonlinear and give rise to significant fre- 
quency shift and mode coupling. We provide analytic ex- 
pressions for the frequency shift, pointing out the crucial role 
played by the anisotropy of the confining potential. For spe- 
cial values of the anisotropy parameter the mode coupling 
is particularly strong and the frequency shift becomes large, 
revealing a peculiar behavior of the Bose-condensed gas. Con- 
sequences on the theory of collapse and revival of collective 
excitations are also discussed. 

3.75.Fi, 67.40.Db 



I. INTRODUCTION 

One of the most important features of an interacting 
quantum many-body system is its response to external 
oscillating fields. The collective modes, which are ex- 
pected to dominate the low frequency response, represent 
a very effective tool for probing the role of interactions 
and testing theoretical schemes. For this reason, mea- 
surements of collective modes in the trapped gases of al- 
kali atoms were carried out soon after the discovery 
of Bose-Einstein condensation 1^ . The remarkable agree- 
ment between measured frequencies and theoretical pre- 
dictions is one of the first important achievements 
in the investigation of these new systems. It provides 
also a clean a posteriori justification of the mean-field 
scheme, based on the formalism of Gross-Pitaevskii the- 
ory I p^ , which is the starting point of most calculations. 
This theoretical approach is expected to give indeed an 
accurate description of the ground state and the excited 
states of such dilute interacting gases at low temperature. 
The same scheme, when the atom-atom interaction be- 
comes dominant compared with the zero-point quantum 
pressure, reduces to the Thomas-Fermi approximation 
for the ground state and to the equations of nondissipa- 
tive hydrodynamics for the excited states ||,|ll|,|l^. This 
is particularly useful for discussing the relevant physical 
properties by means of analytic or semi- analytic results, 
taking advantage of the simplicity and clarity of the hy- 
drodynamic equations. 



The same mean-field approach, which correctly repro- 
duces the frequency of the normal modes of the trapped 
gas in the linear limit (small oscillations around the 
ground state), is suitable to investigating also the non- 
linear dynamics of the systems p^|20[| . The nonlinearity 
is included in the equations of motion through the mean- 
field, which is proportional to the condensate density. 
Thus, measurable effects of nonlinearity could represent 
further clean signatures of Bose-Einstein condensation. 
Among them, large amplitude oscillations of the conden- 
sate can be easily produced in the trapped gases [|l|,D; 
nonlinear effects are expected to give rise to frequency 
shift, mode coupling, harmonic generation and stochas- 
tic motion. 

The purpose of the present work is to derive sim- 
ple differential equations for large amplitude oscillations 
within the formalism of hydrodynamics. In particular 
we will provide analytic formulae for the frequency shift 
of three collective modes in a generic axially-symmetric 
trap. These formulae allow us to discuss the important 
role played by the anisotropy of the trap. We will show 
that special values of the anisotropy parameter (the ratio 
of the axial and radial frequencies of the trap) can be as- 
sociated with strong nonlinear effects even for oscillations 
of relatively small amplitude. The same analytic results 
have interesting consequences in the theory of collapse 
and revival of the condensate. 

In the next section we introduce the basic formalism. 
Then, in section [II, we discuss how the collective modes 
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can be driven and analysed numerically. In section 
we perform a small amplitudes expansion and derive an- 
alytic solutions for the m = and m — 2 modes and 
their frequency shifts. In section ^ we discuss both the 
numerical and the analytic results for different traps. An 
application to the theory of collapse and revival of the 
the oscillations is given in section VI. The paper will end 
with a short summary. 



II. BASIC FORMALISM 

Let us start with the hydrodynamic equations in 
Thomas- Fermi approximation H: 



d_ 

di 



p + V ■ (vp) = 



(1) 



1 



m—v + \/\^Vext+9P+—j^0. (2) 

Density and velocity are related to the condensate 
wave function ^'(r,i) through p = j'l'p and v — 
fi{2mip)~^{^*V^ - ^PV**). The density is normal- 
ized to the number of particles in the condensate, N — 
J drp{r). The external confining potential has the form 
Vextir) = {^/'^)mY,i^'i'rh where = x,y,z. The trap- 
ping frequencies can depend on time, uji = Wi{t), in the 
presence of an external driving force. Their static val- 
ues, LOQi = ti-'i(O), fix the equilibrium configuration of the 
system. For cylindrically symmetric traps one can define 
the radial frequency uj± = luq^ — ujQy and the asymmetry 
parameter A — u!oz/uj±. The external potential provides 
also the typical length scale of the system in each direc- 
tion, = yJhJJjnLjQi) . Since the system is dilute, the 
atom-atom interaction enters only through the quantity 
g = Airh^a/m^ where a is the s-wave scattering length. 

The equations of nondissipative hydrodynamics are 
equivalent to the time dependent Gross-Pitaevskii equa- 
tion for the condensate wave function in the limit 
Na/ano S> 1, if the interaction is repulsive (a > 0). 
In this case, the effects of the zero point kinetic en- 
ergy (quantum pressure) become negligible and the gas 
is dominated by the balance of the internal and external 
potential energies. The stationary solution of the hy- 
drodynamic equations is the Thomas-Fermi ground state 

p^^ir) = \^^^ir)\' =g-'[fi-V,,t{r)] for ^ > K.t(r) 

(3) 

where the chemical potential p is fixed by the normal- 
ization of the density to the number of particle N. The 
hydrodynamic approach works in an excellent way for the 
lowest collective modes of the Sodium atoms trapped at 
MIT where N is of the order of 1 million and more. 
Conversely, it provides only a semi-quantitative descrip- 
tion of the Rubidium gas first trapped at Jila [Q , where 
the number of atoms was smaller (10^-10^). Even in that 
case, however, the measured frequencies converge nicely 
to the hydrodynamic predictions for the largest values of 
N. Compared with the numerical solution of the Gross- 
Pitaevskii equation, the hydrodynamic formalism has the 
advantage of providing analytic results for the dispersion 
law of the collective modes and for other useful quanti- 
ties. 

As already discussed in our previous paper exact 
solutions of the hydrodynamic equations can be found in 
the form 

p{r, t) = ax{t)x'^ + ay{t)y'^ + a,{t)z^ + ao{t) (4) 
V = ^'^[ax{t)x^ + ay{t)y^ + az{t)z^] , (5) 

restricted to the region where p > 0. With this choice, 
equations M-0) transform into a set of coupled differential 



equations for the time-dependent coefficients aj{t) and 
aj{t). One of them is fixed by the conservation of N: 

flo = -{15N/8TT)'^^^{axayazy^^- The equations for the 
others can be further simplified by introducing the new 
variables bi defined by at = —muj^^{2gbxhybzb^)'^. The 
hydrodynamic equations then yield — bt/bi and 

bi + Lofbi ~ ujqJ (bibxbybz) = , (6) 

with i — x,y, z. These equations describe the time evo- 
lution of the widths of the atomic cloud, since the new 
variables bi are directly related to the mean square radii 
and velocities of the system [|l^: bf oc (rf) and bf oc {vf). 
Different derivations of equations (|^) and some applica- 
tions are given in Refs. ||lq-[l7[|. A variational approach 
including the zero point quantum pressure, beyond the 
Thomas-Fermi approximation, has been also presented 
in Ref. |Q. Note again that the frequencies uii, entering 
the second term of (||), can depend on time and, hence, 
these equations can be used for describing time varying 
traps, as well as the expansion of the gas after a sudden 
switching-off of the confining potential. The set of so- 
lutions defined by the scaling transformations (|[^) does 
not exhaust all possible motions of the trapped gas. For 
instance, the motion of the center of mass can be included 
by adding terms linear in x, y or z and other solutions can 
be found including terms of the form xy, xz or yz. How- 
ever, they are well suited to study the collective modes 
of lowest multipolarity and energy, namely the m = 
and m = 2 modes, where m is the azimuthal angular 
momentum in the cylindrically symmetric trap. In the 
following we will apply equation (|^) to these modes in 
the nonlinear regime. 

III. OSCILLATIONS OF A DRIVEN 
CONDENSATE 

At equilibrium one has bi = 1 and bi — 0. One can per- 
turb the system by modulating the trap frequencies for a 
certain time and then let it oscillate freely. Formally, this 
means that equations ^ have to be solved using a time 
dependent frequency of the form ujf{t) = Wgjl -I- 2fi{t)], 
where fi(t) is an appropriate sinusoidal function within 
a finite time interval, while the factor 2 is introduced for 
convenience. An example is given in Fig. 1 for the case of 
a trap with asymmetry parameter A ~ uJoz/^± = %/8, as 
in Jila experiments. Starting from equilibrium, an oscil- 
lation is induced by choosing fx(t) — rjsinQj^t, fy = — fx 
and fz = 0, with a driving frequency fid = V2 uj±. The 
driving force is switched-off at t ~ 20ujJ^. The three 
curves correspond to the width bx{t) plotted for different 
values of the driving strength, 77. The free oscillations are 
undamped, since the theory is restricted to zero temper- 
ature and does not include any dissipation. In the small 
amplitude limit, rj —^ 0, the oscillation coincides with the 
m = 2 normal mode predicted by the linearized hydrody- 
namic equations . It corresponds to a quadrupole-type 
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excitation in the xy-plane. For larger amplitudes, the re- 
sponse of the system is slightly shifted in frequency, as 
can be seen from the figure, and the oscillations are no 
more purely sinusoidal. This is associated with the oc- 
currence of harmonic generation and mode coupling. In a 
similar way, one can excite the to = modes. The one at 
low energy corresponds to an in-phase oscillation of the 
width along x and y and out-of-phase along z. The one at 
high energy is an in-phase compressional mode along all 
directions (breathing mode). In all cases, one can solve 
numerically the equations (O) and extract amplitude and 
frequency of the excited modes. A simple way consists in 
doing a best-fit with a sinusoidal function to the appro- 
priate width bi{t), after the switching-off of the driving 
force. One can also perform a Fourier analysis of the 
signal. We used both methods, finding the same results 
provided the time interval contains at least three or four 
complete oscillations. 




time 

Evolution in time of the radial 



width b: 
1 



FIG. 1. 

a trap with A = \/8. Time is in units of w^^. The os- 
cillation is driven by an external force for the first 20 time 
units; then the condensate oscillates freely. The three curves 
correspond to increasing the strength of the driving force 
(ry — 0.01,0.04,0.06, see text). In the small amplitude limit, 
the oscillation coincides with the m = 2 normal mode. 



decoupled in the linear limit, that is, when the perturba- 
tions ^-,^+ and ^2 are much smaller than 1. Thus, these 
functions represent the normal modes of the system: ^2 
is the m = 2 mode, having frequency \/2 uj^^ while ^_ 
and ^-1- are the low-lying and high-lying m = modes, 
respectively, whose frequencies are ^Jq^ (|] . We note 
also that 7^ for any value of A. 

In Fig. 2 we plot the functions [t] for the oscillation 
of largest amplitude already shown in Fig. 1. The driv- 
ing force is tuned on the to = 2 mode, ^2j but, due to 
nonlinear coupling, the low-lying m = mode, is also 
significantly excited, while ^+ remains very small. The 
coupling between these modes causes the irregular oscil- 
lation of bxj given by (0) and plotted in Fig. 1. These 
oscillations have a very large amplitude (about 80% of 
the radial width); for smaller oscillations the coupling 
with ^_ tends to vanish. For the smallest oscillation in 
Fig. 1 both ^_ and ^+ are negligible with respect to the 
driven mode ^2 and one has also (t) ~ 1 -I- ^2 (i) • 
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FIG. 2. Evolution in time of the three functions ^2 (solid 
line), ^_ (dashed line) and ^+ (dot-dashed Une) in a trap 
with A — \/8- Time is in units of lo^^ . In the limit of small 
amplitude, these functions coincide with the m = 2, low-lying 
m = and high-lying m = normal modes, respectively. The 
driving force is the same as for the largest oscillation in Fig. 1. 



In order to point out the effects of mode-coupling, it 
is convenient to introduce suitable combinations of the 
widths hi{t) as follows: 

6,(t) = i + e-(t)+e+(t)+6(0 (7) 
6,(t) = i + e-(t)+e+(t)-6(0 (8) 

K (i) - 1 + {q- - 4) (t) + {q+ - 4) C+ (i) , (9) 
with 

qzp = 2 + (3/2) ip (l/2)v/9A4 - IGA^ -|- 16 . (10) 

Inserting these combinations into (|^) one obtains equa- 
tions for ^_ , ^+ and ^2 ■ Such equations are completely 



A strong enhancement of nonlinear effects can be ob- 
tained not only by increasing the strength of the driving 
force, but also by changing the anisotropy of the trap. 
In fact, for special values of A = uJz/oja^, frequencies of 
different modes, or of their harmonics, can coincide. An 
example is shown in Fig. 3, where we plot the functions 

[t) for an oscillation driven as in the case of the small- 
est oscillation in Fig. 1. The driving force is again tuned 
on ^2 with a relatively small strength, but now the oscil- 
lation exhibits a strong nonlinear behavior. The crucial 
parameter is the value of A, which in Fig. 3 is chosen to 
be a/16/7. As we will discuss in the next sections, this 
special value of A is associated with a resonance of the 
two modes £,2{t) and S,+ {t)] indeed, the figure shows an 
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evident beating of these two modes. 
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FIG. 3^ Same as in Fig. 2 but for a trap having 

A = ^ 16/7. The beating is between the ^2 mode, which 
is excited by the external force, and the 5+ mode, which is 
coupled through the second harmonic of ^2. The driving force 
is the same as for the smallest oscillation in Fig. 1. 



At this point, before presenting further numerical re- 
sults obtained from equations we prefer to explore 
the first nonlinear corrections to the uncoupled normal 
modes. As we will see, an expansion of at small am- 
plitudes provides useful analytic formulae, which contain 
the main physical ideas and give accurate results. 



The two m = modes are not excited directly, but they 
can be driven by the m = 2 mode via the last term in 
(|T^). One can easily find conditions for resonances: the 
frequency of the m = mode, ,Jq^ wj^, should be equal 
to the frequency of the second harmonic of the m = 2 
mode, 2^/2 lo±. This never occurs for the lowest ra — Q 
mode, since is always smaller than 10/3. Conversely, 
for the high- lying mode a resonance is obtained for 
A ~ -\/l6/7. For this value of the asymmetry param- 
eter, a small driving force, appropriate for exciting the 
m = 2 mode, give rise to a coupled motion of both the 
m = 2 and high-lying m = modes. This is exactly 
what we have already shown in Fig. 3 by solving the hy- 
drodynamic equations for the same special value of 
A. This enhancement of nonlinearity, due to resonances 
via second harmonics, has been receiitly suggested also 
in Ref. Even though equations are valid in 

the limit of small amplitudes, they give the correct con- 
ditions for the occurrence of strong nonlinear effects in 
the exact hydrodynamic solutions. 

The above equations can be also used to predict the 
frequency shift of the collective modes at the lowest order 
in the amplitude A. In fact, equations ( [l^ ) can be solved 
explicitly, by inserting ^2 at the leading order: ^2(i) = 
A2 cos(-\/2 u)±t). One finds 



±{q± - 4) 



2(g+-g_)(g^-8) 



Using these solutions into (|ll|), one gets an equation for 
^2 at the next order: 



IV. ANALYTIC EXPANSIONS AT SMALL 
AMPLITUDE 

Equations (^), with Ui = ujoi and i — x,y,z, describe 
free oscillations in the trap. They can be expanded by 
assuming bi{t) = 1 + 5bi(t) and keeping the lowest orders 
in the small perturbation 6bi. The linear terms provide 
the hydrodynamic spectrum of normal modes, while har- 
monic generation is included in the higher order terms. 
One can also use definitions in order to get three 

coupled equations for the functions ^_ , ^+ and ^2 , up to 
quadratic or cubic terms. We show here an example of 
an explicit calculation for the m = 2 mode, while for the 
m = modes we will present only the final results. 

Let us assume the m = 2 mode be excited in the cylin- 
drically symmetric trap by some suitable driving force, 
in such a way that ,^2(i) ~ v4cos(-\/2 u)±t). If A <C 1, the 
coupling with the other modes is small and, hence, the 
m — oscillations ^_ and ^+ are expected to be of order 

. Then, by expanding up to third order in A, 

one gets the following equations: 

C2 + ^i6[2 - i-q- - (+q+ + 2^2] 







(11) 



2 I q± 



q± \q+ 



(12) 



6 + 2c.i(l-A^i^4^U2 



2ujj 



2 16 - 7A2 
16_3A2 
16 - 7A2 



(14) 



This is a special case of the more general class of equa- 
tions for anharmonic oscillators: 



(15) 



One can easily prove that the solution S^{t) has a fre- 
quency to which depends on the amplitude A through 



UJ = [ 1 



12^4 + SQ2 



By comparing ( [l4| ) and (|5|) one identifies 



n^.^2.Ui-Al'^-^ 



73 = 2col 1 



16_3A2 

16 - 7A2 



(16) 



(17) 



(18) 



and 72 — 0. Inserting these definitions into (|T^), one gets 
finally 
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\/2 



(16 -5A^) 2 
4(16 -7A2) 2 



(19) 



which is the analytic expression of the frequency shift as a 
function of the anisotropy of the trap. Here the pathology 
at A = is evident: as already said, it corresponds 

to a resonance induced by harmonic generation. 
In a more general form, one can write 



w = t^o[l + ^(A)A2 



(20) 



where loq is the frequency of each normal mode in the 
linear regime (v^wj, and ^Jq^ui^ for the m = 2 and 
m = modes, respectively) while A is its amplitude. 
The coefficient for the m — 2 mode is then 



^2 (A) 



(16 -SA') 
4(16- 7A2) 



(21) 



Expansion (^0|) is expected to be reliable whenever 
1 51^2 <C 1 and this excludes values of A too close to 
resonances. The coefficients <5ip(A) for the m = modes 
can be calculated straightforwardly, with the same pro- 
cedure used above for m ~ 2. In particular, one has 
to start exciting one of the two modes, ^ip, in the form 
AzpCOs(y/g^ uj±t); the other one is excited to the order 
A?^, while ^2 is never excited, due to the axial symmetry 
of both the equations and the initial conditions. Expand- 
ing (^ up to terms in A'^ , one gets two coupled equations 
for of the type (p^, while ^2 is decoupled. One finally 
obtains 



6, 



5.2(9±-2)(gT-4)((7^-5) 
2^ 



15 



1 



16 (<7± - 

9 (g± - 4) 
3 g=F - 3 



q±)iq± - q^f 

2 [-q^ + 2\^q^ 



15 A2 

4 q\ 

9A2+8]' 



20 q^ - q± 



[-lOA^g^ + 37A2 + llq^ - 54] . (22) 



Two resonances are found by exciting the low-lying m = 
mode. They occur when the frequency of the high- lying 
mode is equal to the second harmonic of the low-lying 
mode; this happens for A = (\/125 ± \/2^)j\fl2 (i.e., 
A « 0.683 and A « 1.952). Conversely, no resonances 
are found by exciting the high-lying m = mode. As 
an example of numerical values of 5zp, let us consider a 
spherical trap, A = 1; in this case one has q_ = 2, g+ = 5, 
(5_ = 11/12 and 5+ = —1/6. The latter value was already 
used in Ref. |^ to give an estimate of the collapse time 
for the high lying m = mode. 

The main results of these analytic expansions, together 
with the numerical solutions of (|^), will be discussed in 
the following section. 
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FIG. 4. Coefficient of the quadratic expansion ( |2o[ ) for the 
frequency shift of the low-lying m = mode, as a function 
of the anisotropy parameter A = tJoz/ajx- Divergences are 
found for A = (^25 ± 729)/ ^72. The values of A for the 
experimental traps of Refs. are also indicated. 
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FIG. 5. Coefficient of the quadratic expansion ( |20[ ) for the 
frequency shift of the high-lying m = mode. For A = 
the coefficient is about 31. 
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V. FREQUENCY SHIFT AND MODE COUPLING 
FOR DIFFERENT TRAPS 

So far have developed two methods for investigating 
the behavior of the m = and m = 2 modes in a cyhn- 
drically symmetric trap. One is the numerical solution 
of the coupled differential equations ^ for the widths 
bj{t), which corresponds to solve exactly the equations 
of hydrodynamics for these collective modes. The other 
one is an expansion at small amplitude, which provides 
analytic solutions as well as simple formulae for the fre- 
quency shift. 

In Figs. 4-6 we first show the three curves for the coeffi- 
cients ^-(A), (5+ (A) and (52(A), from (^l|-p2|), entering the 
quadratic expansion (|2^) of the frequency shift. These 
curves can be considered the main results of the present 
work. The occurrence of resonances is a peculiar feature 
of the dynamics of the condensate. Close to these special 
values of A where (5(A) diverges, the frequency shift is ex- 
pected to show an anomalous enhancement and the time 
evolution of the widths of the condensate should become 
very irregular, due to significant mode coupling. 
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FIG. 6. Coefficients of the quadratic expansion ( |20[ ) for 
the frequ ency s hift of the m = 2 mode. A divergence is found 
at A = y/l6/7. 



In order to test the accuracy of the small amplitude 
expansion, yielding the curves in Figs. 4-6, we compare 
the results obtained from both the numerical solution of 
and the quadratic law (^) for specific values of A. 
In Fig. 7 we show the results for the m — 2 and low- 
lying m — modes at A = VS, as in Jila experiments 
101 . In this case one finds the coefficients (52(-\/8) = 3/20 
and S-{^/8) = 1.636 x 10^'^. Solid lines are obtained 
from the numerical solutions of (^, using definitions (^ 



g), and performing sinusoidal fits and/or Fourier analysis 
of ^2 and The frequency shift of the m = mode is 
smaller than the one of m = 2 especially for small am- 
plitudes; in the A limit, the agreement between the 
numerical results and the quadratic law (dashed lines) is 
very good. In the m = 2 case it remains good even for 
large amplitudes, while higher order corrections seem to 
be important for the m — mode at amplitudes larger 
than 0.2 (corresponding to relative amplitudes of the or- 
der of 20% in the radial width of the condensate) . In the 
latter case, one notices also that the coefficient 5- is very 
small, since the value A = \/8 is incidentally very close to 
a root of the function (5_(A), as shown in Fig. 4. Thus the 
first nonzero contributions to the shift come from powers 
or higher. We note also that the abscissa in Fig. 7 
is the amplitude of ^2 {t) and (t) and not of the width 
bx{t), which is the observed quantity at Jila. The rea- 
son is that the quadratic expansion ( pO| ) have analytic 
coefficients only for the normal modes However, since 
A = \/8 is relatively far from resonances, the coupling 
between different modes in (^js]) is rather weak and the 
relative amplitude of the oscillations of bj. practically co- 
incides with the amplitude of each normal mode, within 
the range of Fig. 7. 
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FIG. 7. Frequency shift for the m = 2 and low-lying 
m = modes in the Jila trap as a function of their am- 
plitude. Solid line: from the numerical solution of (|^; 
dashed line: quadratic expansion @ with (52 = 3/20 and 
(5- = 1.636 X 10"^ from (|l|- 

It is not easy to compare these results with the avail- 
able data for the frequency shift in the Jila experiments 
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0. The data show no significant frequency shift for the 
low-lying m = mode and a relatively large shift for the 
TO = 2 mode. The fact that the shift is larger for to = 2 
than for m = is in agreement with our predictions, but a 
quantitative comparison should take into account at least 
two further effects. First, the data have been taken after 
switching off the trapping potential. As discussed in our 
previous work (see also |^^), the expansion of the gas 
causes an amplification of the relative amplitude of the 
observed oscillations. This amplification is different in 
the three directions, depending on the anisotropy of the 
trap and on the excited modes. Second, the assumption 
Na/auo 3> 1, which is at the basis of the hydrodynamic 
approach, is not well satisfied by the samples used so far 
at Jila. In this case, the role of the quantum pressure, be- 
yond the Thomas- Fermi approximation could be relevant 
also for the dynamics of the expansion. 

In Fig. 8 we plot the frequency shift of the low-lying 
TO = mode for the MIT trap H with A = 0.077. Again 
the solid line comes from the evolution of ^-(i), calcu- 
lated numerically from (||), while the dashed line is the 
quadratic approximation (pO|), whose coefficient is now 
5- = 1.481. The agreement between the two curves is 
good. In this trap, the condensate is cigar-shaped and 
one measures the oscillations of the axial width bz{t). 
Looking at equation (^), one notes that the relative am- 
plitude of the oscillations of hz is about four times the 
amplitude of the normal mode In fact, for such a 
small value of A, the quantity behaves like (5/2)A^ 
and, hence, can be neglected in (^); moreover, the cou- 
pling with the high-lying m = mode is weak and thus 
one can also neglect ^+ at the lowest order. This means 
that the shift becomes (5-/16)^2, if is the rel- 

ative amplitude of the axial width hz- The shift in fre- 
quency is thus very small; for instance, when the axial 
width oscillates with a relative amplitude of 40%, the 
amplitude of f_ is about 0.1 and the predicted shift is 
less than 2%. Actually the oscillations in Ref. ||^, which 
do not exhibit any frequency shift, have been imaged af- 
ter the free expansion of the condensate and correspond 
to relative amplitudes of the order of 10%, or less, for 
hz before the expansion. Therefore, the predicted shift is 
practically zero within the accuracy of the available ex- 
periments and the law can not be tested. The situa- 
tion can be greatly improved by the use of nondestructive 
imaging techniques, as in the trap of Ref. Q. 

As concerns nonlinear coupling of to = and m ~ 2 
modes, the values of A in the Jila and MIT traps are not 
of particular interest. The curves in Figs. 4-6 suggest dif- 
ferent choices, namely A = y/16/7 for the m — 2 mode 
and A = (%/T25 ± ^29)/^^ for the low-lying to = 
mode. It is also worth mentioning that the curves in 
Figs. 4-6 come from an expansion up to the third or- 
der in A, as in ([Tl[-p^, and the coupling occurs via the 
second harmonics. When the amplitude is large, one ex- 
pects higher order terms to become significant, providing 
other special values of A for mode coupling. For instance. 



the hydrodynamic equations (||) predict a beating of the 
TO = 2 and the high-lying m = modes for A = y/63/11. 
This beating can be calculated analytically by extending 
equations ( pT[]T^ ) to the next order in A; for that value 
of A, one then finds that the third harmonic of ^2 has 
the same frequency of ^+ . Finally, we remark that other 
values of A can give rise to similar effects, even without 
harmonic generation. This happens when an accidental 
degeneracy occurs in the spectrum of normal modes, in- 
cluding |to| > 2, as recently discussed by Ohberg et al. 
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FIG. 8. Frequency shift for the low-lying m = mode 
in the MIT trap as a function of its amplitude. Solid line: 
from the numerical solution of (|^; dashed line: quadratic 
expansion (0) with (5_ = 1.481 from (S). 



VI. APPLICATION TO THE THEORY OF 
COLLAPSE AND REVIVAL 

The fact that the frequency of the normal modes de- 
pends on their amplitude has interesting consequences 
for the phenomenon of collapse and revival of collective 
excitations. The latter belongs to a class of quantum 
phenomena having no classical analogue. Examples are 
the collapse and revival of coherent quantum states ob- 
served in atomic Rydberg wave packets p6t , in molecular 
vibrations p7[ | and for atoms interacting with an electro- 
magnetic field in a resonant cavity |^ (see |^,|o| and 
references therein, for a theoretical discussion) . The same 
concepts have been recently applied to the trapped con- 
densate (see refs. |25 3^ and references therein). In this 
confined Bose system, the collapse of collective excita- 
tions originates from a dephasing during the oscillation. 
In Ref. p3] the dephasing was associated with the fiuctu- 
ations of the particle occupation number and the collapse 



7 



time was estimated within Hartree approximation. One 
of us has recently developed a theory for the same 
process, but based on the fluctuations in the number of 
quanta of oscillation. The main idea is that a collective 
mode is, in general, a coherent superposition of station- 
ary states of an oscillator and the number of quanta of 
oscillations is the natural quantity for classifying those 
states. The estimated time scale for the collapse turns 
out to depend, as a key ingredient, on the frequency shift 
of the normal modes. We will recall here the main steps 
of the theory and then we will use our new results for 
the frequency shifts in order to provide a quantitative 
estimate of the collapse time. 

As we saw in section the collective variable ^ can be 
considered as solutions of the equations for anharmonic 
oscillators. Let us write the frequency, up to the first non 
linear correction, in the form 



UJ = UJQ + SuJ = LUo{l + kE), 



(23) 



where E is the energy of the oscillations. We assume 
here the condition of weak nonlinearity, \k\E <C 1, even 
though the theory can be generalized to the case of large 
amplitude. In order to provide a quantum description, 
one can use the semiclassical expression hcu ~ (dEn/dn), 
which allows one to to rewrite equation ( p3| ) in the quan- 
tum form LUn = En/h = ujqti + bin? /2. Here the quantity 
n is the number of quanta (n 3> 1) in a given excited 
state of energy Em while the coefficient b is given by 
b — TiloqK. The oscillations in the experiments were 
driven through a sinusoidal force into a coherent state of 
the oscillator. The wave function of a such a stationary 
state can be written in the form 



V' = X! exp(-ia;„t) 



where 



I c„ p= — exp(-n) 



exp 



(n — nY 
2fi 



(24) 



(25) 



where n is the average number of quanta (n 3> 1) and 
E — fiLo^n. Let us now consider the oscillator co-ordinate 
£,{t) and calculate its average over the state (p^ps]). If 
one considers only n ^ n ± 1 transitions , the result is 



cos[(a;o + bn)t\. For small enough 



values of i, one can replace the summation over n by an 
integral and one gets a Gaussian damping of the oscilla- 
tion according to (^) '-^ exp(— n6^t^/2) = exp[— (i/rc)^], 
where 



= {n/2Y" \b\= LUoiEhLOo/2y^' \k \ . (26) 

Since the amplitude of the oscillations is proportional to 
VH, the resulting amplitude dependence of the collapse 
time, Tc, is the same as in the theory of Ref. ||3^. The 
periodicity of {^{t)) gives also the revival period — 
27r / (HloqIkI) as in Ref. |^. The meaning of Tc and is 
also shown schematically in Fig. 9. From the expressions 



for Tc and one easily sees that Tc ~ ^y{^/n) ^ Tj.. 
Note that in quantum mechanics the measurement of an 
oscillator coordinate ^ is, generally speaking, destructive 
because it affects the oscillator momentum p^. In other 
words, it modifies the coefficients c„, thus preventing the 
observation of collapse and revival. In order to measure 
properly (^(i)), one must repeat cycles of observations 
with different replica of the system in different instants 
of time. This is, however, exactly the procedure used in 
the available experiments 



exactly 

mm- 



A 
V 




time (arb. units) 
FIG. 9. Schematic picture of collapse and revival. Both 
time and {Q are in arbitrary units and the oscillations are just 
a pictorial view of the phenomenon. The two time-scales, Tc 
and Tr are indicated. 



According to (|2^) the collapse time Tc decreases when 
n increases. However, one has to keep in mind that 
the theory is restricted to the weak nonlinear regime 
\k\E ^ 1, which corresponds to n ^ 1 / {Tiloq\k\) . Thus, 
the coll apse tim e is subject to the condition Tc 3> Tmm = 
1 / {ujQ ^Jhu)[)\K\) and the collapse can be observed only if 
Tmin is shortcr than the typical time scale for the dissipa- 
tive damping of the oscillations, which is always present 
at finite temperature. 

In order to provide quantitative estimates of the col- 
lapse time for trapped Bose condensed gases, one must 
calculate the energy of oscillations, E, as a function of 
their amplitude. This is what we have done in the pre- 
vious sections. One has only to write the proper rela- 
tion between the nonlinearity coefhcients k, and 5. To 
this aim, it is enough to calculate the energy E in lin- 
ear approximation, as twice the mean kinetic energy 
E = m JdrpQV^, where po(r) is the equilibrium density of 
condensate and v(r) is the velocity field associated with 
each normal mode. Direct integration for m — 2 gives: 



E = -fiNAl 



(27) 



where A2 is the amplitude of the mode and /i is the 
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chemical potential of the system, calculated within the 
Thomas-Fermi approximation. Similarly, for the m = 
modes one finds 

E=^ [2q^ + {q^ - 4) iSq^ - 10)] tiN , (28) 

where the quantities q^{X) are defined in (|l^. The above 
energies can be rewritten in a more compact form as E ~ 
ejiA^, where e(A) is the appropriate coefficient for each 
normal mode in (^-28). Inserting this result in ^ 
comparing with the quadratic law (^), one finds 



and 



and revival of the oscillations. By a proper choice of 
the anisotropy parameter A one can in fact increase sig- 
nificantly the coefficient 5 and, consequently, lower the 
collapse time to an observable scale. Obviously, for real 
traps, one should also take into account possible non- 
linear effects originating from trivial nonharmonic cor- 
rections to the magnetic confining potential. But the 
present scenario for collapse and revival seems plausi- 
ble and, in any case, the first step for its confirmation 
should be the observation of the predicted amplitude de- 
pendence of the frequency shift. 



(29) 



The final result for the collapse time, from (pq), is then 



V2 e^iN 



2e 



5 I A I 5 



hujo 



(30) 



Let us estimate for the m = 2 mode, assuming typ- 
ical experimental parameters of Ref. Q]: N = 4500, 
Lu±/2Tr = 132 Hz, A = V8 and a = 110 Bohr radii. 
With these parameters one finds y^JlNjhUZ ~ 195. The 
frequency of the m — 2 mode in the linear limit is 
ujo = V2^±, while the coefficient of the quadratic ex- 
pansion (|20|) is 62 = 3/20. Equation (||) gives €2 = 4/7. 
With a reasonable choice for the amplitude, A = 0.2, the 
final result for the collapse time is Tc = 4.98 s. The au- 
thors of Refs. [p]j3l|] reported a lifetime of the m — 2 mode 
of the order of 100 ms. It means that, under the experi- 
mental conditions, the dissipative damping is too strong 
and the collapse can not be observed. It is not hopeless, 
however, to discover this effect under appropriate condi- 
tions. The point is that, according both to measurements 
pit and to theoretical considerations , the dissipative 
damping decreases rapidly by lowering the temperature 
of the gas. We believe that experiments at lower temper- 
ature would permit to observe the quantum collapse of 
collective modes in these macroscopic objects. It is worth 
mentioning, in this connection, the recent suggestion to 
cool adiabatically the gas by changing the trapping fre- 
quency |34|. 

The low-lying m — Q mode, which also has been 
observed, exhibits a very small frequency shift and is 
not, for this reason, a proper object to search for col- 
lapse, at least within the range of amplitude where the 
quadratic expansion (20) holds. Conversely, the high- 
lying m = mode looks more promising. For this mode, 
using again the parameters of the Jila trap, one has 
LOo ~ 4.98wx, e « 198 and 5+ « 30. The oscillations are 
strongly anisotropic and the amplitude in the z-direction 
is larger that in the radial one. According to (^ one 
gets Az ~ (q+ — 4)A+ « 20.8A+. By assuming again 
Az = 0.2, the collapse time becomes Tc = 1.4 s. 

From these results, it appears evident that the reso- 
nances described in the previous sections of this work 
are very promising even for the observation of collapse 



VII. CONCLUSIONS 

In this paper we have investigated the behavior of col- 
lective excitations of trapped Bose gases in the limit of 
zero temperature. The formalism of hydrodynamic equa- 
tions is suitable for deriving both numerical and analytic 
results, valid in the limit Na/ano ^1- As already 
shown by other authors, those equations describe the 
lowest collective excitations in the linear regime (normal 
modes), in agreement with available experiments. They 
can be also used for the dynamics of the system in nonlin- 
ear regime. We have shown here that they provide non- 
trivial predictions for the frequency shift of the normal 
modes. We have discussed the case of to = and to = 2 
modes in a generic cylindrically symmetric trap and we 
have explored the peculiar behavior of the condensate for 
special choices of A = ujqz/^a^^ the asymmetry parame- 
ter of the trap. For those special traps, the frequency of 
an oscillation becomes equal to the second harmonic of 
another one. This degeneracy enhances nonlinear effects, 
producing significant consequences in measurable quanti- 
ties: for instance, the time evolution of condensate shape 
can show irregular patterns and the frequency shift can 
become rather large. The physical picture has been sup- 
ported by simple analytic results. We have also discussed 
an application to the the theory of collapse and revival 
of the collective oscillations. 
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